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Abstract 

We present a quenched lattice calculation of the lowest order (0(a 2 )) hadronic contribution to 
the anomalous magnetic moment of the muon which arises from the hadronic vacuum polarization. 
A general method is presented for computing entirely in Euclidean space, obviating the need for 
the usual dispersive treatment which relies on experimental data for e + e~ annihilation to hadrons. 
While the result is not yet of comparable precision to those state-of-the-art calculations, systematic 
improvement of the quenched lattice computation to this level is straightforward and well within 
the reach of present computers. Including the effects of dynamical quarks is conceptually trivial, 
the computer resources required are not. 

PACS numbers: 12.38. Gc, 13.40.Em, 14.60.Ef, 14.65.Bt 
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The magnetic moment of the muon is defined by the q 2 — > (static) limit of the vertex 
function which describes the interaction of the electrically charged muon with the photon, 

r,(p2,Pi) = IpFM 2 ) - j ! —{l P i-il P )F 2 {q 2 ), (1) 

where is the muon mass, q = p 2 —pi is the photon momentum, and pi, p 2 are the incoming 
and outgoing momentum of the muon. Lorentz invariance and current conservation have 
been used in obtaining Eq. ^ Form factors F\{q 2 ) and F2(q 2 ) contain all information about 
the muon's interaction with the electromagnetic field. In particular, i*i(0) = 1 is the electric 
charge of the muon in units of e, and g = 2i<\(0) + 2F 2 (0) = 2 + 2 F 2 (0) is the Lande g- 
factor, proportional to the magnetic moment. The anomaly, defined as half of the difference 
of g from its tree level value, which the Dirac equation predicts to be 2 for an elementary 
spin 1/2 particle, is a M = -F 2 (0). Thus, F 2 (0) = at tree level, and corrections to ^(0), 
and therefore a M , start at 0(a) in QED where a = e 2 /4ir is the fine structure constant. 
iq(0) = 1 to all orders due to charge conservation. 

The most precise measurement ever of the muon's anomalous magnetic moment was re- 
cently carried out at Brookhaven National Laboratory In the authors quote three 
standard deviation discrepancies between the Standard Model and experiment . The the- 
oretical and experimental uncertainties are roughly the same and were added in quadrature. 
The dominant theoretical uncertainty resides in hadronic loop corrections arising from the 
hadronic vacuum polarization (0(a 2 )) (see Figure ^) and hadronic light-by- light scattering 
(0(a 3 )), and it is clearly of interest to reduce these errors. Presently, the 0(a 2 ) hadronic 
contribution is calculated by using a dispersion relation and the experimental value of the 
total cross-section for e + e~ annihilation to hadrons to relate the imaginary part of the vac- 
uum polarization to the real part. This calculation is very precise, though a discrepancy with 
a calculation that uses r decay data may indicate a theory error as large as five percent Q] 
and reduces the disagreement with experiment to roughly 1.6 standard deviations. A purely 
theoretical, first principles, calculation has been lacking and is desirable, and also has sev- 
eral advantages over the conventional approach. For instance, the separation of QED effects 
from hadronic corrections is automatic, as is the treatment of isospin corrections if different 
quark masses are used in the simulation. Thus it is possible that lattice calculations may 
eventually help to settle the above mentioned discrepancy between e + e~ annihilation and 
t decay. 
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The method described here is simple and direct. We begin with Ref. Q which describes 
the computation of multi-loop graphs in perturbation theory through the expansion of the 
integrand in terms of hyperspherical polynomials. The key is that the entire integral, includ- 
ing external momenta, can be Wick-rotated into Euclidean space and the angular integrals 
done so that what is left is an integral over the magnitude of the loop momentum. If the 
graph can be set up in a certain way, then after the external momenta are analytically 
continued on-shell, the integral over the loop momenta is performed without distorting the 
integration contour: K 2 = — > oo, K 2 the squared Euclidean loop momentum. For the case 
at hand, this means the photon propagator can only depend on the loop momentum. This 
is important because if the photon propagator depends only on K 2 , then so will the vacuum 
polarization. Thus, the renormalized hadronic vacuum polarization function, calculated on 
the lattice and in Euclidean space, can be directly inserted into the one-loop diagram for the 
anomalous magnetic moment without analytically continuing the vacuum polarization back 
to Minkowski space, or requiring its value in a region not accessible to lattice calculations, 
namely K 2 < 0. The condition on the photon propagator is easily met in this case as can 
be seen from the assignment of momenta in Figure H 

Following (4,!^, to extract F 2 (q 2 ) from Eq. Q apply a projection operator P p and compute 
Tr (P p T p ). Omitting the quark loop for the moment, applying the Feynman rules, and taking 
the limit q 2 — > 0, the diagram in Figure HI gives [2(3] 

,(D_„ a,- r d * k 1 1 A6(p-fc) 2 , 4, 2 



a M = e * / TTrn 77 — 2 — -^71 — - l 2 + o k + 4 (p • k ) > ( 2 ) 

M J (27r) 4 ((p - k) 2 + m 2 - it) 2 k 2 - it \ 3ml 3 / 

where pi^ = p q/2. Analytically continuing the entire amplitude to Euclidean space, 
including the external momenta, performing the angular integrations, and then analytically 
continuing p 2 — > — back on-shell, we are left with 

f°° mlK 2 Z 3 (l- K 2 Z) a f°° 
/ dK Vt \ ran = dK f(^)> 3 

Jo 1 + m l K Z TT Jo 

where Z = —{K 2 — (K 4 + Am 2 l K 2 ) 1 ^ 2 )/2m 2 l K 2 . One can easily verify that this gives the 
leading Schwinger contribution a]u = a/2ir. Because the quark loop does not affect the 
rest of the integral, and more importantly only depends on K 2 , we can simply insert its 
contribution into Eq. El to obtain the 0(a 2 ) hadronic contribution. 

a (2)had = (^_y H dK 2f( K 2)tl(K 2 ), (4) 



a« = - 



where fl(K 2 ) = Air 2 ^\ Q 2 (Jl^K 2 ) — IL,(0)). Qi is the electric charge in units of e and 
Ui(K 2 ) is the vacuum polarization for the i th quark flavor which has been renormalized by 
subtracting its value at K 2 = (from now on we drop the subscript since we work with 
degenerate quarks). In general U(K 2 ) is an analytic function of K 2 . In the conventional 
approach, a]L is calculated via a dispersion relation which relates dt(H(K 2 )) to ^s(H(K 2 )), 
and thus the e + e~ — > hadrons total cross-section. In contrast, here we deal with ^R,(U(K 2 )) in 
the Euclidean, or space-like, region with real K 2 > 0. A simple check of this procedure, which 
is successful, is to insert the one-loop QED vacuum po^fl, andean, — 
to Euclidean space,into Eq. H]and compare the result to the known value of the 0(a 2 ) QED 
contribution to a M |2| - Finally, we note that the kernel f(K 2 ) diverges as K 2 -> 0, so the 
integral in Eq. H] is dominated by the low momentum region. We now turn to the lattice 
calculation of U(K 2 ). 

The vacuum polarization tensor for a single quark flavor with unit charge is defined as 

n^^/dW.^TWfe)) (5) 

where J M (x) = ip{x) r )^ip{x) is the electromagnetic current and () signifies an average over 
gauge and fermion fields. 

In the lattice regularization using domain wall fermionsjsl, 0], current conservation is 
given by A M J M (x) = where A M is the backward difference operator and 

J»{x) = -^2^j(x + fi,s)U\x)(l + Y)^{x,s) -^{x,s)U{x)(l- ^)tp(x + fi,s) (6) 



is the conserved vector current jlO] 2l|. The sum in Eq. |U] is over an extra fictitious 5 



th 



dimension which gives rise to Ad Dirac fermions. These Dirac fermions are chirally symmetric 
up to violations which are exponentially small in the size of this dimension. The Ward- 
Takahashi identity for the two-point function yields 

A"J»{x) (J"(y))t = (7) 

- \ E b[x ~ y) + s ) u \y)( l - ^Wiv* s ) + s)u(y)(i + ^Mv + o, *)) 

s 

+ 8{x-y-v) (${y + z>, s)U\y)(l - j u )i[>(y, s) + ${y, s)U{y){\ + j v )^(y + z>, s)) 

which is valid for each gauge field configuration. The contact terms do not cancel each other 
because J M is a point-split current. After subtracting [ill] 

^ E \ ( + °> s ) U Hv)(l - YMv, s) + ${y, s)U(y)(l + YMv + M ) (8) 
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from the two-point function to cancel the contact terms in Eq. [71 Fourier transforma- 
tion yields the Ward-Takahashi identity g^TP^g) = 0. The polarization tensor calcu- 
lated in Euclidean space, on the lattice, is rP^(g) = (q^q u — q 2 5^ u ) n(g 2 ) which follows 
from Euclidean and gauge invariance. The lattice four-momentum is q = 2/ a sin 
(g M = 27rn fl /N fl , n M = 0, 1, . . . , A 7 ^ — 1). The Ward-Takahashi identity, which is valid before 
averaging over gauge fields, has been verified to numerical precision in the present calcula- 
tion (this is a good check that the computer simulation is correct). After subtracting the 
contribution of unphysical heavy fermions from a 5 th dimension of size L s sites [10], fl{K 2 ) 
differs from the continuum vacuum polarization by terms that are 0(a 2 ). 

We have calculated U(K 2 ) in the quenched approximation using the DBW2 gauge 



action 



12j and valence domain wall fermions. Two values of the gauge coupling were chosen 



which correspond to inverse lattice spacing's (set from the p meson mass), a 1 = 1.3 and 
1.96 GeV (see Q). The lattice volumes studied were A" 3 x N t = 8 3 x 24 (1.3 GeV only) and 
16 3 x 32 (1.3 and 1.96 GeV), corresponding to spatial volume V = (1.2 fm) 3 , (2.4 fm) 3 , and 
(1.6 fm) 3 , respectively. For the domain wall fermions we used L s = 8, domain wall height 



M5 = 1.8, and a single 4g? quark mass m/ = 0.04, or roughly 90 and 120 MeV in the MS 
scheme at [A — 1.3 and 1.96 GeV, respectively 14]. 

Shown in Figure El is the vacuum polarization for the large volume at a -1 = 1.3 GeV, the 
most physically interesting one since we are mainly interested in the small q region. The 
agreement with perturbation theory |l5j| is very good, even down to q 2 ~ 0.5 GeV 2 , until 
finally non-perturbative effects begin to be important. For very large values of q lattice 
artifacts dominate, and the lattice and continuum results disagree. The perturbative result 



is evaluated in the MS scheme at /1 = 1/a = 1.3 GeV with a quark mass m„ « 90 MeV 14] 

n 

and has been shifted by a constant [22] in Figure |2 for comparison. We also show the small 
volume results which indicate finite volume effects are negligible until q 2 ^ 0.5 GeV 2 . The 
results from the 1.96 GeV lattice are quite similar, indicating that non-zero-lattice-spacing 
errors are small. 

Since we are interested in the low q 2 region, to use Eq. 0J we fit the lattice data to a simple 
polynomial in q 2 , which allows for a smooth interpolation of the data. This amounts to a 
Taylor expansion about q 2 = 0. Lorentz covariance (really, hyper-cubic symmetry) requires 
Il(g 2 ) depend only on q 2 , and if the quark mass is non-zero, it must be regular as q 2 — > 0. 
A more sophisticated anastz is clearly desirable, especially one that includes finite volume 
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effects. However, the construction of such an ansatz, for example in chiral perturbation 
theory, is outside the scope of this work. The value of n(g 2 ) between q 2 = and the smallest 
value calculated on the lattice is extrapolated from the fit, so it is important to have values 
of q 2 that are close to zero. The smallest value in this study, which is on the large volume, is 
<7min ~ 0.065. Specifically, we use a four parameter fit function, H(x) = a +ai x+a 2 x 2 +a 3 x 3 
where x = q 2 , and fit the data in the range x < 2 GeV 2 . The fit is uncorrelated but performed 
under a jackknife procedure. The final results are insensitive to the number of parameters 
used and the fit range, so long as the chosen combination accurately reproduces the data. 
That is, the larger the range, the more parameters needed to accurately represent the low 
momentum region (see Figure |2J). 

The fit function is next plugged into Eq. 0] and integrated numerically, using MATH- 
EMATICA, up to some cut, q 2 ut . The perturbative value of U(K 2 ) is then used in the 
integral from q 2 ut to oo, after it has been shifted by a constant to match onto the subtracted 
lattice result. The final result is quite insensitive to the value of the cut since f(K 2 ) in 
Eq. 0] is sharply peaked at zero, and with the statistical error on the present data such 
as it is, the perturbative contribution can be ignored entirely since it adds, depending on 
the cut, (1 — 10) x 1CT 10 to aj^ had . However, as the lattice results become more accurate, 
these contributions will have to be carefully included. To be specific, results are given for 
q 2 ut = 1.5 GeV 2 , though a much smaller value yields essentially the same answer. 

Using the above procedure and including degenerate u, d, and s quarks we find a^ ha " d = 
460(78) x 10~ 10 on the large volume (the error is statistical). This is roughly 2/3 of the 
value computed using the dispersive approach[2], and given the approximations in this first 
calculation: quenching, finite volume, and unphysically large quark masses, quite encour- 
aging. We note with respect to the quenching systematic error, the above result is quite 
reasonable: 72% of the dispersive result comes directly from the p resonance pj, essentially 
all that is included in the quenched case. The result depends heavily on the low q 2 region, so 
the final statistical error is still rather large since only a small number of configurations were 
used (27) to compute averages over the gauge-field. The smaller volume result (a -1 = 1.3 
GeV, 138 configurations) is roughly = 318(69) x 1CT 10 , indicating large finite vol- 

ume effects. The 1.96 GeV lattice, which corresponds to a somewhat larger volume, gives 
fl (2)had _ gyg^gg^ x iq-io configurations), in between the large and small 1.3 GeV lattices. 

Since the hadronic loop is given by a sum over all possible values of the momentum of 
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the quarks, one should worry that lattice artifacts due to finite L s may be noticable. Such 
effects were observed in domain wall fermion calculations of the chiral condensate and weak 
16 . \v\ and are well understood. In the latter case, the physical matrix 



matrix elements 



element was obtained by taking a slope with respect to quark mass which easily removed 
the lattice artifact. A similar mechanism is at work here. II (q 2 ) depends on L s through a 
constant, momentum independent shift which vanishes as L s — > oo, but the renormalized 
fl(g 2 ) is independent of L s to a very high degree. This was checked by calculating a]L ^ with 
L s = 4 on the small volume, 1.3 GeV lattice. While the values of 11(0) differ by roughly 30%, 
a£ ? d = 317(64) x 10~ 10 calculated with L s = 4 is in excellent agreement with the L s = 8 
value. Being purely an issue of the physics near the ultraviolet cut-off where the heavy 5d 
fermions play a role (apart from the small additive quark mass, m res h|, which affects the 
low energy physics), this could have been anticipated from Figure 121 Since the agreement 
with perturbation theory is so good, visible effects would only show up as a constant. 

We note that the vacuum polarization can also be fit to a form given by the operator 
product expansion though being a short distance expansion, it is valid for q 2 ^> 1. In 
particular, in this (quenched) study we find no evidence for power corrections beyond the 
operator product expansion, which is in agreement with [ill ]. 

We are optimistic about the prospects for improving this calculation. Further quenched 
calculations on larger volumes will pin down finite volume effects, reduce the statistical 
error (there is room for significant improvement), and probe the chiral limit for the light 
quarks. The disconnected diagram that does not vanish in the non-degenerate case will 
also be computed. This should reduce the error on the quenched calculation below five 

n 

percent. Recent simulations with dynamical domain wall fermions |18| may be valuable in 
estimating the quenching error. Perhaps even more interesting in the near term, recent 
2+1 flavor dynamical fermion calculations using improved Kogut-Susskind fermions I9I ] 
offer a promising avenue to address all the main systematic errors in this first calculation. 
The procedure set down here goes through in just the same way, but with dynamical gauge 
fields replacing the quenched ones. It is the generation of gauge fields with the dynamical 
fermion effects, especially for light quarks, that is so costly. Though we are still a ways 
from competing with the precision of the dispersive method, lattice calculations should 
eventually rival that precision. We have already mentioned that the lattice method avoids 
the problems of disentangling QED effects from the hadronic corrections, and isopin breaking 



7 



can be handled simply by using different quark masses for the u and d quarks. In any case, 
a completely theoretical, first principles calculation, though challenging, is well worth the 
effort. 
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FIG. 1: The lowest order hadronic contribution to the muon anomalous magnetic moment a . The 
muon has outgoing momentum p + q/2 after scattering from a photon with momentum q. The loop 
momentum is k. The blob represents the non-perturbative hadronic vacuum polarization. 

"There is an additional diagram consisting of two quark loops, each with a single insertion of the electro- 
magnetic current, connected by three gluons at lowest order in QCD. This diagram vanishes when summed 
over three degenerate quarks u, d, and s with Q = 2/3, -1/3, and -1/3, respectively which is the case studied 
here. Therefore, its contribution is expected to be small, though this will be checked in future calculations. 
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FIG. 2: The unsubtracted vacuum polarization for the o^ 1 = 1.3 GeV lattice. The two volumes 
described in the text are shown, 8 3 (squares) and 16 3 (circles). The solid line denotes a three-loop 
perturbation theory calculation |l£j], which has been shifted by a constant for comparison, and the 
dashed line the fit of the larger volume data points. 
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